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Unifying Geometric Entanglement and Geometric Phase in a Quantum Phase 

Transition 
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Geometric measure of entanglement and geometric phase have recently been used to analyze 
quantum phase transition in the XY spin chain. We unify these two approaches by showing that the 
geometric entanglement and the geometric phase are respectively the real and imaginary parts of a 
complex- valued geometric entanglement, which can be investigated in typical quantum interferome- 
try experiments. We argue that the singular behavior of the complex-value geometric entanglement 
at a quantum critical point is a characteristic of any quantum phase transition, by showing that the 
underlying mechanism is the occurrence of level crossings associated with the underlying Hamilto- 
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Quantum phase transitions (QPTs) are qualitative 
changes in the properties of many-body systems driven 
by quantum fluctuations at zero temperature. A key 
feature to understand the nature of QPTs is the quan- 
tum correlations between the system degrees of freedom, 
which at the critical point bring about a modification of 
the ground-state level degeneracy. While QPTs have tra- 
ditionally been characterized by appropriate correlation 
functions [HQ, many researchers have recently addressed 
this problem from a quantum information perspective. 
One idea along this line is to use measures of entan- 
glement to characterize QPTs. Indeed, there has been 
a great deal of work analyzing properties of quantum 
entanglement in many-body systems undergoing QPT 
0, Q . Another idea is based on the fact that non-trivial 
geometric phases (CPs) are strongly affected by quan- 
tum and classical correlations residing in the many-body 
quantum states [j| , and especially by the level degenera- 
cies associated with QPTs in such systems. Therefore 
GPs can also be used to detect the presence of a QPT. 
Relations between GP and QPTs have been established 
theoretically and experimentally in NMR 



In this Letter, we examine the geometric measure of 
entanglement (GE) 0, E3] and GPs in the vicinity of 
QPT in a XY spin-chain. We determine a relation be- 
tween these geometric objects, allowing us to identify 
them as the real and imaginary part of a single mea- 
surable quantity, which we call the complex- valued GE. 
We demonstrate that the complex GE is accessible in 
interferometry experiments. The established relation is 
valid for a general quantum many-body system. At a 
quantum critical point both real and imaginary parts of 
the complex-valued GE display the same singular behav- 
ior, which in turns is closely associated with the singular 
behavior of the quantum geometric tensor triggered by a 
level degeneracy. This provides a universal approach to 



the study of quantum critical phenomena. 

Consider the one-dimensional XY model system with 
N sites in a transverse magnetic field. The corresponding 
Hamiltonian reads 
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with periodic boundary condition such that the first and 
(N + l)th sites arc identified. Here, r is the anisotropy 
parameter, h is the magnetic field strength, and a^, 
k = x,y,z, are the standard Pauli operators of the jth 
spin. In the thermodynamic limit (N —> oo), the system 
described by Eq. ([1]) undergoes a QPT at h = he = 1. 
The full phase diagram can be found in Refs. [H, Ell- 
in the anisotropy range < r < 1 the system be- 
longs to the Ising universality class. The isotropic case 
r = corresponds to the XX model, which belongs to the 
Berezinsky-Kostcrlitz-Thouless universality class. The 
standard procedure to solve the eigenvalue problem of 
H(r, h) is to convert the spin operators into fermionic op- 
erators, using successively Jordan- Wigncr, Fourier, and 
Bogoliubov transformations [12J . The ground state in the 
even fermion number sector can be expressed in terms of 
fermionic modes c, as 



\ip(r,h)) = [ cos 6> m (r, h) 

m=0 

+i8mMr,fc)44-m-i]lt-t> (2) 



with tan 29. 
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j(r, h) = r sin 
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n i v" ^ n > ra- 
r < 1 the ground-state in the thermodynamic 

limit is doubly degenerate in the ferromagnetic regime 

h < 1 and singly degenerate for h > 1, with the critical 

value for h c = 1 being a point of conical intersection. 
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Entanglement of the ground state \i/j(t, h)) can be mea- 
sured by approximating it by the closest pure product 
state UJM- This state is given by the maximal overlap 



A, 



Ar, h) = max |($|^(r, h))\ 
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over all pure product states $. The resulting num- 
ber A max (r, h) is the entanglement eigenvalue of the XY 
ground state [iJJ]. GE of the ground state is quantified 
via 



Ei og2 [ip(r, h)] = - log 2 A 2 max (r, h), 



(4) 



which in the thermodynamic limit is characterized by the 
geometric entanglement (GE) density 



s(r,h) = lim £n{t, h), 



(5) 



where £N{r,h) = E\ og2 {ip(r,h)]/N is the entanglement 
per site. Using the translational symmetry of \tp{r, ft,)), 
the closest pure product state for each value of r and h 
takes the form [J] 
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where the parameter £ is chosen so as to maximize 
the overlap, i.e., A max (r,h) = max e \($(£)\i/j(r, h))\ = 
\(^{S,niax)\^(r, h))\. By using Eq. (O, one obtains the 
GE density [J] 



s(r,h) 



max 
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d[i ln[cos 9{[i, r, h) cos 2 (£/2) 



+ sin 9(fi, r, h) sin 2 (£/2) cot tt/j], 
with tan2#(/i, r,h) = r sin27r/i/ (h — cos 27r/i). 
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FIG. 1. (Color online) Derivative of the ground-state GE 
density e(r, h) of the XY model as function of the magnetic 
field h in the thermodynamic limit. Solid curve: Ising limit 
(r = 1); dashed curve: anisotropic (r = 0.5) XY model; dot- 
dashed curve: anisotropic (r = 0.05) XY model. The inset 
corresponds to the Ising limit for different finite lattice size 
N = 11, 18,31,64. 

It has been shown Q that features of quantum critical- 
ity, such as universality, critical exponents, and scaling 
are captured by the singular behavior of the GE of the 
XY ground state. The critical exponents for different 
universality classes have been found by scaling analysis 



of divergences at the singular points of the GE density. 
The singular behavior of the GE density at the critical 
point h c = 1 is shown in Fig. [TJ where (?£jv(r, h)/dh has 
been plotted as function of magnetic field h for different 
values of r in the thermodynamic limit. In the inset we 
plot the same quantity for Ising model for finite values of 
N. The cusp at hf > 1 for < r < 1, coalescing with the 
singularity at h c — 1 as r — > 0, is an interesting feature 
related to properties of the closest product state <5(£ ma x)- 
For each < r < 1 there is a finite value of h = hf for 
which £ max vanishes and f(( max ) is polarized in the z- 
direction 0] . This state remains the closest pure product 
state for all h > hf , causing a discontinuity jump in the 
second derivative of e(r, h) at h = hf. 

Recently, a generic connection between GPs and QPTs 
has been identified 0] . A scaling analysis of GP of the XY 
ground state reveals similar information about quantum 
criticality as the GE density. Top panel in Fig. [2] shows 
the behavior of the GP per site fi 9 N {r, h) of the XY ground 
state, accumulated by adiabatically rotating each spin 
around the z axis via 
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where U tot {<j)) = ®V =1 Uj(4>) 
ent lattice size N. /3 9 (r, h) = lim^v-^oo @%(r, h) is the GP 
density, i.e., the GP per site in the thermodynamic limit. 

The lower panel of Fig. [5] displays the behavior of the 
difference between the GPs per site associated with Cg T ' h:> 
and the evolution 
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U tot (m(r,h))(<l>(r,h)\uLM (9) 
of the closest pure product state |$(r, h)) = |$(£max))- 
Cp ,h ^ is a path of pure product states along which each 
instant state is the closest pure product state of the cor- 

(r h) 

responding instant ground state along the path Cg . 
The entanglement is fixed along this pair of paths. The 
associated GP difference is in this sense related to the en- 
tanglement residing in the initial ground state \tp(r,h)). 

In the thermodynamic limit, we obtain the GP densi- 
ties 

P a {r,h)= lim p%(r,h) 
r-i/2 



= 2n sin 2 8([i, r, h)d[j,, 
Jo 



P(T,h) = lim p%(r,h) = 27rsin 2 (£ max /2) (10) 

N— >oo 

corresponding to Cg 7 "'^ and Cp h \ respectively. 



The GP density difference is 

A/3(r, h) = P 9 {r,h) - P p {r,h). 



(11) 



For < r < 1, Af3(r,h) displays the same (Ising- 
universality class) critical behavior at h = h c as e(r,h). 
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FIG. 2. (Color online) Top panel: Derivative of the GP den- 
sity, GP per particle, corresponding to cyclic evolution of the 
XY ground state given in Eq. {SJ, as function of magnetic 
field h in the thermodynamic limit. Lower panel: h deriva- 
tive of difference between GP densities corresponding to the 
evolutions given in Eqs. © and ©, plotted as function of 
magnetic field h in the thermodynamic limit. Solid curve: 
Ising limit (r = 1); dashed curve: anisotropic (r = 0.5) XY 
model; dot-dashed: anisotropic (r = 0.05) XY model. The 
insets correspond to Ising limit for different finite lattice size 
N = 11,18,31,64. 



Interestingly, even the cusp in e(r, h) at h > h c is cap- 
tured by A/3(r, h) (now in the form of a discontinuous 
jump in the first derivative of A/3(r, /i)), while the ground 
state GP density f3 9 (r, h) is smooth on this h interval. 

For the XX model (r = 0), we observe that A/3(0, h) = 
and therefore A/3(0, ft,) = for the specific evolution 
operator Utot{<t>) UM- O n tne other hand, as indicated 
in the lower panel of Fig. ® ^-A/3(r, h) is singular at 
h = h c = 1 when r^0. Thus, 



lim 



<9A/3(r, h) 



dh 



dA/3(0,/i) 



h=h c 



dh 
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h=h a 



which shows that the XX criticality is detected by non- 
smooth contractible phase difference. Arguments along 
this line have been used in Refs. 0, EH to locate the 
criticality of XX model by using GP. From the analysis 
of Ref. [4J, we see that also de(r, h)/dr scales differently 
for r ^ and r = near the critical point. Just as in the 
case of GE density, the nature of Eq. (|T2"T) identifies the 
different universality classes of XY spin chain. 

We now give a unifying operational interpretation of 
GP and GE for the XY system, in the context of interfer- 
ometry. Let |0), |1) span the state space of an auxiliary 
qubit, playing the role of the two interferometer arms. 
Prepare the initial state = ^(|0) + \l))\tp(r,h)). 
In the |1) arm, project onto the product state $(£)) 
given in Eq. The two spin chain states \ip(r, h)) 



and |$(£)) are exposed to the parallel transporting uni- 
tary operators £/J(0) = ®^ ie I^K- W^k*!^./*))) 
and Ul{<p) = ®^ ie ^K-<*(€)kjl*(€)» ; rcS p C ctively, 

where 4> e [0, n]. Note that U'l((f>) and u\(<p) is the same 
unitary up to an overall phase factor, which assures that 
\(^(^)\U^U^\^(r, h))\ for given r and h is a function of 
£ only. The resulting state U^(tt)\^(^)) in the |1) arm 
is parallel transported along the shortest geodesic back 
to |$(£)), while a U(l) shift e~ if is applied to the |0) 
arm. Finally, the two arms are brought back to interfere, 
resulting in the final state 



A<P a -f) 



(13) 



with ip 9 and <p p being the GPs associated with the evo- 
lution of \ip(r, h)) and $(£)), respectively. The intensity 
in the |0) arm becomes 
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where 



A(r,h;0 = {iP(r,h)\e-^ P m))m)\e^ H'(r,h)) 



\mmr,h))\ 2 e^. 



(15) 



Here, \A(r,h;£)\ = \{*(£)Mr, h))\ 2 and ar gJ 4(r, ft; £) = 
Aip characterize the visibility and phase shift of the in- 
terference fringes obtained by varying /. Note that 
A(r, h; £) is unchanged under any local gauge transfor- 
mation and reparametrization of the paths. The entan- 
glement eigenvalue A max and the GP difference A<p(r, h) 
associated with Cg rJl ^ and C^'^ can be read-out from 
the interference fringes by tuning the single parameter £ 
to the value £ max giving the maximal visibility. 

Following the procedure leading to GE, we introduce 
a complex- valued GE being defined as the following ex- 
tensive entanglement sensitive quantity 

£fo g2 W, h)} = - log 2 A(r, h- £ max ) 

= Eu ^WM-i^£ L - (16) 

The properties of the quantum critical point is character- 
ized by both the real and imaginary parts of the complex- 
valued GE density 



N 



e{r,h) 



; A/3(r,ft) 
1 In 2 ' 



(17) 



The GE and GP densities are thus two different sides of 
the same coin: the complex- valued GE density s c (r,h). 
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Let us now generalize the above ideas to a generic N- 
body system. Assume that the system is prepared in the 
state 14") and thereafter evolves around a loop Cm gener- 
ated by a one-parameter family of local unitary operators 
U(t) = ®^Uj{t) with Uj(t) acting on subsystem j. The 
closest product state |$) evolves under the same U(t). 
Note that the resulting C$ is in general an open path. 
Let Atpi^) = tp(^) - be the GP difference asso- 

ciated with paths Cm and Cm- By following the above 
interferometer scheme, the interference fringes are char- 
acterized by 



A(*) = A^(*)e iA ^) 



(18) 



with |<E>) chosen to maximize the visibility. Similar to 
the XY-case discussed above, A(^) can be used to define 
complex-valued GE of the iV-body system prepared in 
\f r . Furthermore, Afif) can be interpreted in terms of 
the quantum geometric tensor T(v,u) [l5|. This tensor is 
defined asT(v,u) = (v\(l— \ip)(ip\)\u) for tangent vectors 
u and v at a state \4>)- We can write T = j + iF/2, where 
7 is the quantum metric tensor providing the Fubini- 
Study measure of distance in state space and F is the 
symplcctic curvature 2-form responsible for the GP. We 
have 



1 

l -2 



7, 



Acp(V) = f F, 
Js 



(19) 



where G is a geodesic in state space connecting ^ and 
and S is the oriented surface with boundary dS = 
C^*G*G^ 1 *C^ 1 *G- 1 . The inverse denotes the reverse 
direction along the corresponding path. This procedure 
is illustrated in Fig. |3l 




FIG. 3. Paths Cm and C<s> of many-body state $ and its 
closest product state respectively. G is a geodesic in state 
space projecting *]/ onto $. 

Now, assume that the above procedure is applied to 
the family of ground states |*) = |i/>o(A)) of the Y-body 
Hamiltonian H (A) for each parameter value A G M.. One 
can see that on the ground state manifold the quantum 
geometric tensor takes the following form [l6[ 

T _ (^o(A)|^g|^(A))(^(A)|^g|V> (A)) 



[E (X) - E n (X)]< 



where ip n (X) are eigenstates of the system and the in- 
dices /i, v — 1, . . . , dimA'f are labeling the coordinates of 



M.. Note that the imaginary part of is the Berry 
curvature on M. 17|. Since QPTs occurs at level cross- 
ing or avoided level crossing, the energy denominators 
E n (X) — E (X) in Eq. ([20)) show that in the thermody- 
namic limit, 7 = ReT and F = 2ImT will show a singular 
behavior at the QPT. This singularity is reflected in A max 
and Aiy5 through the expression given in Eq. (|19[) . There- 
fore it can be captured by both GE and GP difference, 
which are respectively the real and imaginary parts of 
the complex-valued GE. So the connection between GE, 
GPs and QPT is in fact a generic one. 

In summary, we have established a unifying connection 
between the GP and the GE in a generic many-body sys- 
tem. They are respectively the real and imaginary part 
of a generalized complex-valued geometric entanglement 
and they both become singular at a QPT. We have used 
the quantum criticality of the XY model to derive this 
property, but shown that this result holds in general for 
any QPT. The approach presented here can be tested 
experimentally in an interferometry setup, where the 
geometric entanglement would then yield the visibility 
of the interference fringes, whereas the geometric phase 
would describe the phase shifts. The resulting interfer- 
ence fringes would therefore display large fluctuations at 
the critical point. Test-bed experiment to demonstrate 
the relevance of the complex-valued geometric entangle- 
ment in translational symmetric few-qubit ground states 
would be within reach in, e.g., NMR systems. 
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